POWER TRADING RISK MANAGEMENT SYSTEM 



CORRESPONDING PATENT APPLICATION 

This application is based upon and claims the 
5 benefit of priority from the prior Japanese Patent 
Applications No. 2003-024676 which was filed on January 
31, 2003. 



BACKGROUND OF THE INVENTION 
10 The present invention relates to a power trading 

risk management system for managing a risk in a power 
trading . 

An object of the present invention is to provide a 
method for allowing an electric power company possessing 
15 numerous power plants to procure electric power from a 
spot market on economically favorable terms and to respond 
to numerous electric power demands through related 
financial product transactions on economically favorable 
terms . 

20 Another object of the present invention is to 

provide a power trading risk management technique which 
can manage a risk in a power trading by treating electric 
power in a spot market as a portfolio. 

In Japan, power trading has been liberalized for 

25 retail exclusively for large-volume customers (users of 
20000 V or 2000 kW and above) since March 2000. This 
partial liberalization is to be reviewed after three years 
and the liberalization is expected to be significantly 
expanded in 2004 or thereabout. When power trading is 

30 liberalized as described above, electricity prices will 
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fluctuate depending on market situations. Fluctuations of 
electricity prices cause risks of fluctuations in profits 
and costs to electric power companies, electric power 
brokers and customers. Such situations are similar to that 
5 in the present stock market where stock prices fluctuate. 
In this concern, the financial and securities industries 
have developed methods for hedging (reducing or 
eliminating) risks of price fluctuations by use of 
financial engineering. Such financial engineering methods 

10 are deemed to be also effective in risk hedging against 
electricity prices. 

Financial engineering provides the hypotheses that a 
future price fluctuates entirely at random and that the 
future price is not predictable at all except for a drift 

15 term which corresponds to an interest rate. However, the 
electricity prices are influenced by supply-demand 
relations more significantly in comparison with stocks and 
the like. Meanwhile, since the electricity demands 
generally drop at midnights or on weekends, fluctuations 

20 in a daily cycle or a weekly cycle is observed in the 
electricity prices. Moreover, as a result of large 
fluctuations of demands attributable to weather, the 
prices also tend to fluctuate in accordance with weather. 
Accordingly, seasonal fluctuations in a yearly cycle are 

25 also observed. Since such fluctuations are predictable to 
some extent, the hypothesis that the future price is not 
predictable at all is not always true unlike the case of 
stocks. Meanwhile, when the supply-demand relations are 
strained, the prices may soar and show spike-like 

30 fluctuations. These relations between the demands and 
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prices described above constitute a concept unexpected in 
financial engineering. In addition, the electric power as 
a commodity has various characteristics including that it 
is difficult to preserve, that it generally costs very 
5 expensive, and that it is necessary to use a power 
transmission system to deliver the commodity. Therefore, 
the supply- demand relations are defined by electrical 
rules. From the reasons described above, the fluctuations 
of electricity prices show a slightly different aspect 

10 from the fluctuations of stocks. 

It is necessary to evaluate numerically a quantity 
of risk in order to manage the risk of power trading. To 
this end, it is necessary to model fluctuations of a 
future electricity price. In such a case, the geometrical 

15 Brownian motion model is generally used in financial 
engineering. Here, a conventional technique will be 
explained using the case of a stock as an example. In 
financial engineering, a small deviation in stock price dS 
is generally described as in the following equation (1.) : 

dS 

— = fiat + adz (i) 
2 0 S 

Here, S is a stock price, \x is a drift rate (a trend 
term) , t is time, a is volatility, and z is a variable 
following the Wiener process. 

The volatility is a factor showing uncertainties of 
25 future price fluctuations. This volatility is used in 
financial engineering to describe the magnitude of market 
price. The volatility corresponds to a standard deviation 
calculated on a yearly basis and is defined as the 
following equations (2) : 
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a = s 1 4r 




M , = Info /.?,._, ) J 

Here, Si is a stock price (an electricity price in 
the case of electricity) at time i , and ui is a continuous 
compound interest (or a rate of return) between time i-1 
5 and i. If the price represents a price on every other 
day, then Ui is equivalent to a daily rate of return. 
Meanwhile, if the unit of the time period r is given in 
years, the a becomes volatility of a yearly rate. The 
volatility is a factor indicating the scale of the price 
10 fluctuations . 

In the meantime, the Wiener process is one of the 
Markov's stochastic processes, and is used in the field of 
physics to express a motion of micro particle, or the 
Brownian motion. The Wiener process is expressed as the 
15 following equation (3) : 

dz = £ 4dt ( 3 ) 

Here, dz is a change of z during an infinitesimal 
time period dt , and e is a random sample from a standard 
normal distribution (at an average of 0 and a standard 

2 0 deviation of 1) . The change dz is independent of 
different infinitesimal time periods dt . 

The Wiener process (the Brownian motion) , which has 
a drift term and in which a coefficient of dz is not 1 as 
shown in the equation (1) , is so-called the generalized 

25 Wiener process (or the Ito process) . The equation (1) 
also shows that a logarithm of the stock price performs 
the Brownian motion. Such a stochastic process is called 
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the geometrical Brownian motion. In other words, the 
equation (1) is used to approximate fluctuations in the 
logarithms of the stock prices by using the sum of the 
random fluctuation term that fluctuates according to a 
5 normal distribution. This reflects the fact that 

investors are interested in rates of return more than the 
prices as they are. In evaluation of a risk of stock 
assets, usually, price fluctuations are modeled by the 
above-described geometrical Brownian motion model, and 
10 risk evaluation (measurement) is carried out according to 
the price distribution obtained as a result of such 
modeling . 

For electricity prices, however, there has been a 
difficulty in accurately expressing the price distribution 

15 by use of the geometrical Brownian motion model. 
Hereinbelow, drawbacks of the conventional technique will 
be described based on examples of electricity prices from 
market data at California Power Exchange (CalPX) in the 
United States and at Leipzig Power Exchange (LPX) in 

2 0 Germany. Note that the electricity price data used herein 
are disclosed on the Internet (www.ucei.berkeley.edu/ 
ucei) . 

FIG. 1A and FIG. IB show electricity prices in daily 
average values in the day-ahead market at California Power 

25 Exchange (CalPX) in 1999 and 2000, respectively. 
Meanwhile, FIG. 1C shows electricity prices in daily 
average values at Leipzig Power Exchange (LPX) in 2001. 
FIG. ID shows a transition of closing stock prices of 
Company A in 2001, which is provided herein for reference. 

30 . The first thing apparent from these graphs is that 
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the fluctuations of the electricity prices and daily rates 
of return thereof are considerably greater than the 
fluctuations of the stock prices. As for the volatility, 
in contrast to 55% in the stock prices of Company A, the 
5 CalPX shows 343% in 1999 and 456% in 2000; meanwhile, the 
LPX shows 588% in 2001. Therefore, in each case, the 
volatility of the electricity prices is greater than the 
volatility of the stock prices by nearly 1 digit. 

Prices of many derivative financial instruments are 

10 sensitive to volatility values. Accordingly, the 

conventional financial engineering method has a difficulty 
in hedging a risk of electricity trading. For example, a 
problem of an excessively high option price has been 
pointed out (Yoshiki Murakami et al . : Drafts for Summer 

15 Conference 2002, The Japanese Association of Financial 
Econometrics and Engineering, 2002) . 

Secondly, there is observed a phenomenon that a 
distribution of the rates of return deviates largely from 
the normal distribution. FIG. 2A shows a distribution of 

20 the rates of return in the daily average electricity 
prices at the CalPX, and FIG. 2B shows a distribution of 
the rates of return on a daily basis of the stock prices 
for Company A. Dot lines in the respective drawings' show 
the normal distribution. In the probability distribution 

25 of the ' rates of return at the CalPX in FIG. 2A, it is 
apparent that bottom portions are greater than that of the 
normal distribution. Such a phenomenon is called a fat 
tail. This fat tail is observed frequently in energy- 
related commodities. Although the fat tail is also 

3 0 observed in the stocks, this phenomenon tends to be more 
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significant in the electricity prices. 

A deviation of the probability distribution from the 
normal distribution is expressed by stochastic quantities 
called skewness a 3 and kurtosis a 4 . These quantities are 
5 defined by the following equations (4) : 



Here, x is an average value of data Xi, N is the 

number of the data, and a x is a standard deviation of the 

data. The kurtosis is equal to 3 in the case of the 

10 normal distribution, and becomes greater when the 



0 in the case of a symmetrical distribution. 

In the case of the CalPX electricity prices, the 
skewness a 3 is equal to 0.26 and the kurtosis a 4 is equal 

15 to 6 . 2 . In other words, the distribution is substantially 
symmetrical but is more acute than the normal 
distribution. In the case of the stock prices of Company 
A in FIG. 2B, the skewness a 3 is equal to 0.2 and the 
kurtosis a 4 is equal to 3.3, which is substantially 

2 0 approximate to the normal distribution. A conventional 
financial engineering model usually assumes the normal 
distribution as the distribution of the rates of return. 
Accordingly, the conventional financial engineering method 
cannot be applied directly to power trading. 

25 There is a method of numerically evaluating a 

volatility risk of a market price, which is called a value 
at risk (VaR) . FIG. 3 is a graph showing a loss amount in 
a case where a value of an electricity asset is decreased 




distribution becomes more acute. The skewness is equal to 
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to X L i or below at a certain probability on the assumption 
that the value of the electricity asset fluctuates in 
accordance with the normal distribution based on the 
average jj. and the standard deviation a. This loss amount 
5 is a quantity of risk which is generally called the VaR. 
Here, the electricity asset is equivalent to a product of 
the electricity price and an amount of electric energy. 
It is possible to apply a similar definition in a case of 
multiple assets. Moreover, in a case of different kinds 

10 of assets, it is also possible to apply a similar 
definition by converting the assets into prices. As for 
the multiple assets, if the assets are correlated to one 
another, a standard deviation of a distribution of the 
whole assets can be calculated by use of correlation 

15 coefficients. In FIG. 3, on the assumption that the value 
of the asset is decreased to X Li by 1% probability, the 
loss M-X L i is equal to 2.33xcr according to the normal 
distribution. As for the standard deviation a, a standard 
deviation for a probability distribution in a targeted 

20 future period (in a month, for example) is used therein. 
This value is in proportion with a square root of the time 
in case of assuming the Brownian motion. Accordingly, the 
value can be calculated in an equation "a x = (volatility of 
yearly rate) x (period expressed in yearly unit) 0 ' 5 ". 

25 The case in FIG. 3 is expressed as u the VaR in a 

month is /i-X L i yen at 99% confidence level" . The magnitude 
of the risk is evaluated by the magnitude of this value. 

Even if the distribution is not the normal 
distribution, a similar definition is applicable. FIG. 4 

30 illustrates an example in which an asset value fluctuates 
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in accordance with a probability distribution having a fat 
tail. A point corresponding to a cumulative probability 
(an area ratio of a portion of a probability distribution 
function below X L2 ) of 1% is found and that point is 
5 defined as X L2 . In this case, since X L2 <X L i / the loss 
amount or the VaR becomes greater at the same confidence 
level of 99%. In other words, evaluation of the VaR in 
the distribution having the fat tail on the assumption of 
the normal distribution results in underestimation. For 
10 this reason, it is important for risk evaluation to obtain 
the probability distribution of the future price as 
precisely as possible. However, the conventional method 
has a difficulty in evaluating the probability 
distribution . 

15 In addition, an error in terms of the loss amount is 

another problem at a practical level. FIG. 5A and FIG. 5B 
collectively show a difference between the evaluation 
according. to the VaR and actual profits. This approach is 
called a back test. In this case, daily, losses or gains 

20 ■ (differences between current-day prices and preceding-day 
prices) calculated based on the CalPX electricity prices 
in 1999 and the VaR figures calculated on the assumption 
of the normal distribution are illustrated. FIG. 5B is a 
partially expanded view of FIG. 5A. In the drawings, the 

25 values indicated with bars are the daily losses and gains 
per megawatt hour. The values indicated with dotted lines 
are the VaR figures per megawatt hour on the following 
days at 95% confidence level. Here, the volatility was 
calculated by use of the data for two weeks (14 days) 

30 before the point of evaluation. 
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From the drawings, it is apparent that a frequency 
of losses in excess of the VaR figures is about 1/20 (5%) . 
However, it is learned that the loss amount may be 
enormous once the loss is incurred. Therefore, although 
5 the conventional method can accurately evaluate the 
frequency of losses, the method would underestimate the 
loss amount . 

In addition to the foregoing problems, the risk 
management is far more complicated in the case of the 
10 electricity prices as compared to the stock prices because 
the relations with the demand or the weather must be taken 
into account . 

Furthermore, the business strategy of an electric 
power company is drastically changed by the liberalization 

15 of power market. Specifically, major uncertainties prior 
to the liberalization were a demand fluctuation and an 
unexpected stoppage of a power generator (fuel costs have 
been basically passed to selling prices) . Therefore, an 
object of the risk management has been to forecast the 

20 demand and to minimize power generation costs. In 
contrast, after the liberalization of power market, major 
risk factors are attributable to the fluctuations of the 
electricity prices. . However, as described later, such 
risks are basically advantageous to the electric power 

25 company when appropriately managed. Specifically, the 
electric power company will have an option "not to 
generate power" after the liberalization. Such an option 
does not necessarily mean abandonment of responsibility of 
power supply, but means an option to stop the power 

3 0 generator and to procure electricity from a market for 
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sales instead. In this case, a proportion between an 
amount of power generation by the electric power company 
and an amount of procurement from the market becomes 
important. A brief example will be described below in 
5 this concern. 

Before the liberalization, assuming that a power 
generation cost per unit electricity is C, a fixed 
electricity price per unit electricity is P, and an amount 
demanded is L, then a profit of the electric power company 

10 per unit electricity is expressed by G 0 = L(P-C). 

Meanwhile, a profit of power generation after the 
liberalization will be as follows. Here, in order to 
discuss an ideal case, an assumption is made that there 
are sufficiently fluid spot markets. In addition, 

15 technical restrictions of the power generator are ignored 
and it is assumed that the power generator can freely 
generate power within the maximum electric power Qmax- 
Based on these assumptions, if the power generation cost 
per unit electricity is C, the fixed electricity price per 

20 unit electricity is" P, the amount demanded is L (L<Qmax) , 
an amount of electricity generated is Q (Q<Qmax) / an amount 
of procurement from the market is B (which represents an 
amount of sales to the market when B is negative) , and a 
spot price per unit electricity is S, the profit of power 

25 generation is described as follows: 

(1) if S<C, the amount of electricity generated Q=0, the 
amount of procurement B=L, and the profit G=L ( P- S ) ;>G 0 ; and 

(2) if S>C, the amount of electricity generated Q=Qmax/ 
the amount of procurement B=- (Qmax-D , and the profit 

30 G=L(P-C) + (Q MAX -L) (S-C)>G 0 - 
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That is to say, it is possible to gain the profit by 
stopping the power generation and procuring the whole 
amount from the market when the spot price is lower than 
the power generation cost. In other cases, it is possible 
5 to gain an additional profit by selling a portion 
exceeding the amount demanded out of the capacity of 
electric power facilities in the market. In each case, it 
is apparent that the profit is increased more after the 
liberalization. In reality, it is necessary to find B for 
10 maximizing the profit while satisfying the technical 
restrictions. Meanwhile, if there is a cost C which 
arises even when the power is not generated, then a 
condition of the procurement from the market is expressed 
by S<C-C . 

15 In addition, when emission trading for C0 2 gas and 

the like is put into practice in the future, there is a 
possibility of an additional increase in profit by means 
of selling an emission right when stopping the power 
generation and procuring the electricity from the market 

20 instead. Since a thermal power plant has such an emission 
right as a matter of course, it is possible to gain profit 
by selling the right. 

In consideration of one power plant and the spot 
market as described above, it is obvious that an optimum 

25 combination between these factors exists. In the case of 
using multiple power generators, a method of deciding the 
optimum combination is far more complicated because the 
power generation costs vary among the respective power 
generators . 

3 0 The electric power company today owns many power 
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plants with different power generation costs, and 
therefore decides the optimum combination of the power 
plants to meet the demand which is determined by external 
factors. However, when the power trading is liberalized, 
5 a method of deciding the optimum combination between the 
power plants and the demand is naturally changed due to 
the reasons that the electricity price varies with time 
and that the price varies depending on each demand. A 
portfolio concept used in the financial engineering field 

10 is applicable to such a method. The portfolio concept is 
a method of combining assets with different risks (time 
variation of the prices) and different profits most 
suitably. For example, it is possible to reduce the risk 
by combining assets of different types. Although a 

15 holding period is not considered very much in terms of 
simple portfolio optimization, there are concepts of a 
period for transmitting the electricity and a period for 
receiving the electricity in the case of an electricity 
asset. Therefore, the portfolio optimization method is 

20 required to take account of the time. 

Various types of the electricity prices exist even 
today depending on conditions of contracts. The degree of 
freedom of price setting is increased when the power 
trading is liberalized. Accordingly, the electric power 

25 company will have many types of electricity assets with 
different prices and different holding periods. Assuming 
that the power generation cost is constant, the 
electricity assets will have different rates of return and 
different maturities. Derivative financial instruments 

30 such as futures or options are deemed to appear for the 
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electricity assets sometime in the future. Methods for 
managing such derivative financial instruments are yet to 
be established even in overseas countries where the 
liberalization is advanced. 
5 One portfolio management method considering the 

holding period is asset and liability management (ALM) . 
When a mismatch between short-term funds and long-term 
funds occurs in a financial institution in terms of both 
fund operation and fund procurement, profitability 

10 declines due to differences in interest rates. To avoid 
this phenomenon, various methods have been developed to 
control the assets and the liabilities comprehensively. 

The history of the ALM started a long time ago in 
early 1980' s, involving maturity ladder analysis (for 

15 grasping repayment trends depending on interest rate 
ranks), maturity gap analysis, and sensitivity analysis to 
begin with. Later, as the importance of present value 
risk management has been recognized, and the ALM has 
developed into duration analysis, the value at risk (VaR) , 

20 earning at risk (EaR) , and the like. Recently, the ALM is 
also applied to credit risk management. The usual ALM 
mainly considers interest-rate fluctuations. However, 
this method is also applicable to the risk management of 
power trading. Nevertheless, the ALM method has not been 

25 used for management of electricity assets to date, except 
for some cases such as the VaR. 

Accordingly, the electricity obtained from the power 
plant and the electricity procured from the market are 
deemed as financial assets with fixed maturities and 

3 0 dividend- periods, and electricity demands based on various 
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contracts are deemed as financial liabilities with fixed 
maturities and dividend periods in consideration of 
contracted amounts of electric energy, contract periods, 
contract conditions, and the like. Now, consideration 
5 will be made on a method of managing the financial assets 
and the financial liabilities as a portfolio while 
considering the maturities and the dividend periods 
thereof . 

When the electric power facilities for supplying the 

10 electricity are managed depending on the magnitude of the 
risks so as to form the portfolio together with the 
electricity demands, the method for the optimum 
combination of the power generators (an economic dispatch) 
also becomes different. A combination for minimizing the 

15 power generation costs has been conventionally selected 
upon decision of the power generators to be operated 
(costs for starting and stopping the power generators are 
ignored here for the purpose of simplification) . Such an 
approach corresponds to decision of a combination of power 

20 generator outputs to minimize the fuel costs with respect 
to the demand encountered. For example, assumptions are 
made that there are two power generators, a required 
amount of power generation is P 0 , and optimum generated 
outputs are Pi and P 2 , respectively. When the power 

25 generation costs of the respective power generators are 
expressed by fi(Pi) and f 2 (P2) / then it is a question of 
deciding the values of Pi and P 2 to minimize f 1 (Pi) +f 2 (P2) 
on the condition of P 0 =Pi + P2- The optimum solutions for 
this question can be found by use of the Lagrange 

30 undetermined multiplier method. Here, with respect to 
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F=f x (PJ +f 2 (P 2 ) +X (P 0 -Pi-P 2 ) , the following 
equations (5) need to be resolved: 



simultaneous 



dP l dP x 

8F 3f 2 



A = 0 



(5) 



Therefore, the following equation (6) is obtained 
ultimately : 

dfx _ ¥2 



(6) 



dP x dP 2 

The X here is called an incremental fuel cost. That 
is, even when there are many power generators, those power 
generators should be appropriately operated so that the 
incremental fuel costs become equal . 

For example, if fi(Pi)=aiP! 2 and f 2 (P2) =a 2 P 2 2 , then the 
solutions of the equation (6) are as follows: 



R = 



F = 



(7) 



Nevertheless, the final goal for the electric power 
company is not to minimize the power generation costs but 
to maximize the profit. When the electricity price is 
constant as in the past, minimization of the power 
generation costs is equal to maximization of the profit. 
However, things are different when the electricity price 
fluctuates as a result of the leberalization of power 
market. The method similar to the above is still 
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applicable when considering only the individual power 
generators even in the case where the electricity price 
fluctuates. Actually, the minimization of the power 
generation costs is not always equal to the optimum 
5 combination of the power generators when considering the 
costs for starting and stopping the power generators and 
the like. In general, however, the costs for starting and 
stopping the power generators are relatively small. 

The above -described simplification is not applicable 

10 when the prices vary depending on the demands. Moreover, 
considering that a power generator with a small risk 
should meet a demand with a small risk, for example, the 
minimization of the power generation costs is not equal to 
the optimum combination of the power generators because 

15 the fluctuations of the targeted electricity prices are 
different. Meanwhile, the minimization of the power 
generation costs is not equal to the optimum combination 
of the power generators when a power generation plan is 
made by selecting the power generators in response to the 

2 0 periods and the amounts of the demands from the market and 
so on. Therefore, the economic operation method based on 
the incremental fuel costs as in the past does not provide 
the optimum result under the liberalization of power 
market . 

2 5 As described above, when power trading is 

liberalized, more various electricity commodities will 
appear than those at present. In the meantime, the 
electricity prices will vary depending on the demands and 
the electricity prices even vary with time. To deal with 

3 0 such fluctuations, the conventional technique has a 
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difficulty in comprehensively managing the assets 
including the power supply and the power demand. 
Moreover, from the viewpoint of the electricity asset, 
time factors such as the electricity transmission period 
5 or the electricity reception period are important. 
However, since the conventional method does not take such 
factors into account, it is difficult to optimize the 
portfolio. In addition, since the price distribution of 
such assets is largely deviated from the normal 

10 distribution, the conventional model causes a large error. 
On the other hand, when the electric power company 
executes power trading under the liberalized electricity 
market, the electric power company cannot maximize the 
profit by merely considering the power generation costs. 

15 Furthermore, it is difficult to maximize the profit while 
maintaining the risk within a tolerance, and it is 
difficult to optimize the risk management for each power 
plant . 

2 0 SUMMARY OF THE INVENTION 

One aspect of the present invention is a power 
trading risk management system, in which electricity 
obtained from a power plant and electricity procured from 
a power trading market are deemed as financial assets with 
25 fixed maturities and dividend periods, while electricity 
demands based on various contracts are deemed as financial 
liabilities with fixed maturities and dividend periods in 
consideration of contracted amounts of electric energy, 
contract periods, contract conditions, and the like, 

3 0 whereby these assets are comprehensively managed as an 
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electricity portfolio considering the maturities and the 
dividend periods . 

Another aspect of the present invention is a power 
trading risk management system comprising: an electricity 
5 procurement planning unit for producing an electricity 
procurement plan by combining electricity to be generated 
by one or plural owned power generators and electricity to 
be procured from a market; a generation-procurement curve 
producing unit for producing a generation-procurement 

10 curve based on the electricity procurement plan; a 
portfolio producing unit for producing a portfolio of 
electricity to be generated by owned power plants and 
electricity to be procured from a power trading market 
that is matching to the generation-procurement curve ; a 

15 risk evaluation unit for evaluating a risk of the 
portfolio; a profit estimating unit for estimating a 
profit of electricity sale according to the portfolio; a 
portfolio reorganizing unit for reorganizing the 
portfolio; a best portfolio proposing unit for judging a 

20 best portfolio which can maximize the profit while 
maintaining the risk in the profit of electricity sale in 
a certain period within a tolerance among the reorganized 
portfolios. 

According to this aspect, the power trading risk 
25 management system can further include: means for deciding 
a combination of power generators to be operated to 
maximize a profit accrued from power generation. 

According to this aspect, the power trading risk 
management system can further include: an estimation of 
30 future electricity demand unit for estimating fluctuations 
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of future electricity demand according to past electricity 
demand fluctuations; and an estimation of electricity 
price fluctuation unit for estimating future electricity 
price fluctuations according to the past electricity 
5 demands, past price fluctuations and a relationship 
between electricity demand and price in a predetermined 
period as well as the estimated fluctuations of future 
electricity demand; and in the system, it is possible to 
configure that the portfolio producing unit includes a 

10 price of an emission right for carbon dioxide in the 
portfolio; and the best portfolio proposing unit judges 
the best portfolio which can maximize the profit while 
maintaining the risk in the profit of electricity sale in 
a certain period within the tolerance among the 

15 reorganized portfolios which are including the price of 
the emission right of carbon dioxide. 

According to this aspect, the power trading risk 
management system can further include: an estimation of 
future electricity demand unit for estimating fluctuations 

20 of future electricity demand according to past electricity 
demand fluctuations; and in the system, it is possible to 
configure that the portfolio producing unit includes a 
financial product related to a weather in a corresponding 
region in the portfolio; and the best portfolio proposing 

25 unit judges the best portfolio which can maximize the 
profit while maintaining the risk in the profit of 
electricity sale in a certain period within the tolerance 
among the reorganized portfolios which are including the 
financial product related to the weather. 

30 According to this aspect, in the power trading risk 
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management system, it is possible to configure that the 
risk evaluating unit manages a position and calculates a 
risk index for the electricity portfolio by use of the 
relationship between demand and price which varies 
5 depending on a country, a region, and time wherein the 
system is operated. 

According to this aspect, in the power trading risk 
management system, it is possible to configure that the 
risk evaluating unit uses at least any of volatility, risk 

10 sensitivity, skewness of a rate of return distribution, 
kurtosis of the rate of return distribution, a percent 
point of the rate of return distribution, a percent point 
of a price distribution, a value at risk and an earning at 
risk for managing the position and calculating the risk 

15 index for the electricity portfolio. 

According to this aspect, in the power trading risk 
management system, it is also possible to configure that 
the risk evaluating unit uses a probability distribution 
different from a normal distribution as a distribution of 

2 0 the rate of return attributable to a power trading upon 

risk evaluation of the portfolio. 

According to this aspect, in the power trading risk 
management system, it is possible to configure that the 
risk evaluating unit uses a probability distribution 
25 different from a normal distribution that is calculated 
from a financial Boltzmann model as a distribution of the 
rate of return attributable to a power trading upon risk 
evaluation of the portfolio. 

Another aspect of the present invention is a power 

3 0 trading risk management system, comprising: a power 
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generation risk parameter evaluating unit for simulating a 
fluctuation of a profit of each power plant and evaluating 
a risk parameter of power generation by use of a 
fluctuation of a fuel price; an electricity procurement 
5 risk parameter evaluating unit for evaluating a risk 
parameter of electricity to be procured from a power 
trading market; an electricity demand risk parameter 
evaluating unit for evaluating a risk parameter of an 
electricity contract with each customer; a portfolio 

10 producing unit for producing a portfolio according to a 
proportion of electricity to be generated by owned power 
plants and electricity to be procured from a power trading 
market, the risk parameter of power generation, the risk 
parameter of electricity procurement and the risk 

15 parameter of electricity contract; a risk value evaluating 
unit for evaluating a risk of the portfolio; a portfolio 
reorganizing unit for reorganizing the portfolio by 
adjusting the proportion of the electricity to be 
generated by owned power plants and electricity to be 

20 procured from the market to maximize a profit while 
maintaining the risk amount within a tolerance; and an 
output ting means for deciding the proportion of the 
electricity to be procured from the market which can 
maximize the profit while maintaining the risk amount 

25 within the tolerance as an optimum combination and 
outputting the optimum combination as a power generation 
plan. 

According to this aspect, the power trading risk 
management system can further include: means for deciding 
30 a combination of power generators to be operated to 
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maximize a profit accrued from power generation. 

According to this aspect, in the power trading risk 
management system, it is possible to configure that the 
risk value evaluating unit uses an asset and liability 
5 management method of any of maturity ladder analysis, term 
gap analysis, and duration gap analysis for management of 
the portfolio. 

Still another aspect of the present invention is a 
power trading risk managing method, comprising the steps 

10 of: producing an electricity procurement plan by combining 
electricity to be generated by one or plural owned power 
generators and electricity to be procured from a market; 
producing a generation-procurement curve based on the 
electricity procurement plan; producing a portfolio of 

15 electricity to be generated by owned power plants and 
electricity to be procured from a power trading market 
that is matching to the generation- procurement curve ; 
evaluating a risk of the portfolio; estimating a profit of 
electricity sale according to the portfolio; reorganizing 

20 the portfolio; and judging a best portfolio which can 
maximize the profit while maintaining the risk in the 
profit of electricity sale in a certain period within a 
tolerance among the reorganized portfolios. 

Still another aspect of the present invention is a 

25 power trading risk managing method, comprising the steps 
of: simulating a fluctuation of a profit of each power 
plant; evaluating a risk parameter of power generation by 
use of a fluctuation of a fuel price; evaluating a risk 
parameter of electricity to be procured from a power 

30 trading market; evaluating a risk parameter of an 
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electricity contract with each customer; producing a 
portfolio according to a proportion of electricity to be 
generated by owned power plants and electricity to be 
procured from a power trading market, the risk parameter 
5 of power generation, the risk parameter of electricity 
procurement and the risk parameter of electricity 
contract; evaluating a risk of the portfolio; reorganizing 
the portfolio by adjusting the proportion of the 
electricity to be generated by owned power plants and 

10 electricity to be procured from the market to maximize a 
profit while maintaining the risk amount within a 
tolerance; deciding the proportion of the electricity to 
be procured from the market which can maximize the profit 
while maintaining the risk amount within the tolerance as 

15 an optimum combination; and output ting the optimum 
combination as a power generation plan. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1A is a graph showing electricity prices in 
2 0 daily average values at California Power Exchange (CalPX) 
in 1999. 

FIG. IB is a graph showing electricity prices in 
daily average values at California Power Exchange (CalPX) 
in 2000. 

25 FIG. 1C is a graph showing electricity prices in 

daily average values at Leipzig Power Exchange (LPX) in 
2001 . 

FIG. ID is a graph showing a transition of closing 
stock prices of Company A in 2001. 
30 FIG. 2A is a distribution graph showing rates of 
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return in the daily average electricity prices at 
California Power Exchange (CalPX) . 

FIG. 2B is a distribution graph showing rates of 
return on a daily basis of the stock prices for Company A. 
5 FIG. 3 is a graph showing a loss amount (VaR) when a 

value of an asset is decreased to X Li or below at 1% 
probability on the assumption that the asset fluctuates in 
accordance with a normal distribution based on an average 
|_i and a standard deviation o. 
10 FIG. 4 is a graph showing a loss amount when a price 

is decreased to X L2 or below at 1% probability on the 
assumption that the price fluctuates in accordance with a 
probability distribution having a fat tail. 

FIG. 5A is a graph showing profits of daily average 
15 electricity prices at California Power Exchange (CalPX) 
and showing VaR figures on the following days at 95% 
confidence level. 

FIG. 5B is a partially expanded view of FIG. 5A. 

FIG. 6 is a block diagram showing a functional 
20 configuration of a first embodiment of the present 
invention . 

FIG. 7A is an explanatory view for a portfolio of 
electric power demands. 

FIG. 7B is an explanatory view for a portfolio of 
25 conventional electricity procurement. 

FIG. 7C is an explanatory view for a portfolio of 
electricity procurement according to the first embodiment. 

FIG. 8 is a table showing a concept of electricity 
asset and liability management. 
3 0 FIG. 9 is a block diagram showing a functional 



configuration of an electricity procurement plan producing 
unit of a second embodiment of the present invention. 

FIG. 10A and FIG. 10B are a flowchart showing a 
process carried out by the electricity procurement plan 
5 producing unit of the second embodiment. 

FIG. 11 is a block diagram showing a functional 
configuration of an electricity procurement plan producing 
unit of a third embodiment of the present invention. 

FIG. 12 is a flowchart showing a process carried out 
10 by the electricity procurement plan producing unit of the 
third embodiment . 

FIG. 13A is a graph showing a relation between an 
electricity demand and an electricity price at California 
Power Exchange (CalPX) in February, 1999. 
15 FIG. 13B is a graph showing a relation between an 

electricity demand and an electricity price at California 
Power Exchange (CalPX) in August, 1999. 

FIG. 14A is a graph showing a relation between an 
electricity demand and an electricity price at Leipzig 
20 Power Exchange (LPX) in February, 2001. 

FIG. 14B is a graph showing a relation between an 
electricity demand and an electricity price at Leipzig 
Power Exchange (LPX) in August, 2 0 01. 

FIG. 15 is a block diagram showing a functional 
25 configuration of a fourth embodiment of the present 
invention . 

FIG. 16 is a flowchart showing a process carried out 
by the electricity procurement plan producing unit of the 
fourth embodiment . 
30 FIG. 17 is a view for removing periodicity out of 
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electricity prices. 

FIG. 18A is a graph plotting variation with time of 
skewness by using data of California Power Exchange 
(CalPX) . 

5 FIG. 18B is a graph plotting variation with time of 

kurtosis by using the data of California Power Exchange 
(CalPX) . 

FIG. 19A is a graph showing fluctuations of 
electricity prices close to a normal distribution, which 
10 are calculated by use of a financial Boltzmann model 
according to a fifth embodiment of the present invention. 

FIG. 19B is a graph showing fluctuations of 
electricity prices largely deviated from the normal 
distribution, which are calculated by use of the financial 
15 Boltzmann model according to the fifth embodiment of the 
present invention . 

FIG . 19C is a graph showing daily rates of return 
corresponding to the fluctuations of electricity prices 
shown in FIG. 19A. 
20 FIG. 19D is a graph showing daily rates of return 

corresponding to the fluctuations of electricity prices 
shown in FIG. 19B. 

FIG. 20A is a graph showing a difference between a 
normal distribution model and the financial Boltzmann 
25 model employing VaR figures at 95% confidence level. 

FIG. 20B is a partially expanded view of FIG. 20A. 

FIG. 21A is a graph showing a relation between an 
absolute value of a daily rate of return and a demand at 
the CalPX. 

30 FIG. 21B is a graph showing a relation between a 
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daily rate of return and a traded volume of stocks for 
Company A. 

FIG. 22 is a graph showing a relation between an 
absolute value of a daily rate of return on a current day 
5 and a daily rate of return on the preceding day at the 
CalPX. 

FIG. 23 is a graph showing a result of fitting by 
use of a temperature function. 

FIG. 24 is a block diagram showing a functional 
10 configuration of a sixth embodiment of the present 
invention . 

FIG. 25 is an explanatory view showing an output 
screen of a comprehensive power trading risk management 
system incorporating various functions of the first to 
15 sixth embodiments of the present invention. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Now, embodiments of the present invention will be 
described in detail with reference to the accompanying 

20 drawings. FIG. 6 shows a functional configuration of a 
power trading risk management system according to a first 
embodiment of the present invention. In this power 
trading risk management system, a demand data inputting 
unit 1001 inputs targeted demand related data, which are 

25 individual amounts demanded as a time function and risk 
indices thereof. A demand curve producing unit 1002 
produces a demand curve 1010 for a certain period by use 
of the demand related data inputted by the demand data 
inputting unit 1001. A total generation-procurement curve 

30 producing unit 1003 produces a total generation- 
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procurement curve 1020 proximately. The demand curve 1010 
to be produced by the demand curve producing unit 1002 is 
equal to the total generation-procurement curve 1020 to ' be 
produced proximately by the total generation-procurement 
5 curve producing unit 1003. 

An electricity procurement plan producing unit 1004 
matches amount and time of power generation (or 
procurement) as closely as possible with amount and time 
of a demand, and thereby produces a plan for operating 

10 power generators and a plan for procuring electricity from 
a market or another electric power company as shown by a 
power generation-procurement plan curve 1030. Moreover, 
the electricity procurement plan producing unit 1004 
outputs the plans as target data 1005 for generation and 

15 procurement. FIG. 7A to FIG. 7C show a result of the 
above -described procedures in comparison with a 
conventional method . 

FIG. 7A to FIG.7C show a principle of producing an 
optimum power generation plan by use of the power trading 

20 risk management system of the first embodiment. In 
general, a maximum amount of electric energy and a 
delivery period are set forth in an electricity contract. 
Accordingly, a discussion in terms of a certain section of 
time is not appropriate herein, but a discussion in 

25 consideration of a time axis is required. FIG. 7A 
conceptually illustrates this approach. In this case, 
there are five customers, and a sum of demands of the five 
customers constitutes the demand curve 1010. The lateral 
axis indicates the time. Actual contract conditions are 

30 not simple; in addition, electricity to be used is not 
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always constant at any time. However, for the purpose of 
simplification, the individual demands are represented in 
rectangular shapes. On the contrary, a procurement side 
conventionally operates power generators in the order of 
5 low power generation costs as shown by a total power 
generation-procurement curve 1040 in FIG. 7B. When a 
power generator having the lowest power generation costs 
reaches the maximum generation capacity, then a power 
generator having the second lowest power generation costs 

10 is operated, and so on. In reality, a base power source 
is sometimes selected due to a difficulty in starting and 
stopping operation or changing output. Accordingly, it is 
not always true that the power source having the lowest 
costs is simply the first to be operated. Nevertheless, 

15 when the price does not fluctuate with time, which means 
when there is no risk, the above-described concept 
generally works out . There is no reason for applying an 
expensive power source deliberately as a base power 
source . 

20 On the contrary, when the electricity price 

fluctuates with time along with the liberalization of 
power market, it is not appropriate to decide the power 
generators to be operated merely depending on the power 
generation costs. In other words, it is necessary to 

25 consider risk factors. As rates of return vary among the 
power generators and price volatility risks are also 
different, it is necessary to adopt a so-called portfolio 
approach. A rate of return Y in a portfolio is equal to a 
weighted average of rates of return yi in terms of 

30 respective assets. Specifically, when a weight of the 
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respective assets is wi, the rate of return Y is given by 
the following equation (8) : 

y = 2>,. x . (8) 

i 

In this case, a variance cr Y 2 of Y is expressed by the 
5 following equation (9) : 

€T r 2 =Y. w i 2<T ? ■^mi w i w j <7 ij (9) 

i i j*i 

where, a±j is a covariance between the assets i and 
j. Assuming that there are i-1 pieces of assets, 
estimation of an aspect of change in risk, which is 

10 variance of risk, in a new portfolio caused by adding one 
new asset i is the most important part of the portfolio 
theory. The aspect of change in risk in the entire 
portfolio to be caused by adding the i-th asset thereto 
can be obtained by partial differentiation of the equation 

15 (9) with respect to w ± . Specifically, the following 
equation (10) is applied hereto: 

<^ = 2w i a i 2 +2Yw j cr iJ =2Cov(y,,y) (10) 

The equation (10) corresponds to a covariance 
between the i-th asset and the entire portfolio. 

20 Specifically, it is apparent that the risk in the new 
portfolio is not influenced by the variance of the i-th 
asset but by the covariance with other assets. In short, 
an electric power company needs to judge as to whether it 
is appropriate to acquire a new demand or not on the basis 

25 of a relation with other demands. 

The above-described value at risk (VaR) is also used 
when evaluating the risk in the portfolio. For example, 
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the VaR of a current portfolio is assumed to be 100 
million yen. Here, an assumption is made that there is a 
necessity to enter into a new electricity contract 
estimating a profit of 20 million yen. The value of this 
5 transaction is assumed to be 30 million yen, for example. 
In this case, a proportion between the risk and the return 
is 3:2, and the transaction is not deemed efficient. 
However, if the entire VaR is increased by 10 million yen 
and consequently reaches 110 million yen by incorporating 

10 the new transaction into the portfolio, then the 
transaction is efficient because the proportion between 
the risk and the return is 1:2. Thus, it is not 
appropriate to judge a benefit of a transaction merely by 
considering the profits of the respective power trading. 

15 The same holds true with a relation between a supply 

and demand. As the electricity contract is also a sort of 
properties, the contract is divided into an asset, which 
is a positive property, and a liability, which is a 
negative property. When applying this concept to a 

20 practice in the financial and securities industries, the 
asset corresponds to the electricity to be generated or 
the electricity to be procured from the market, and the 
liability corresponds to the electricity to be supplied to 
the customers. 

25 Usually, the asset is equal to a sum of the 

liability and a capital. However, the capital is assumed 
to be 0 in this case. Imbalance between the demand and 
supply in the power trading leads to a blackout. 
Accordingly, the asset and the liability are deemed to be 

3 0 balanced without any margin that corresponds to the 
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capital. An ALM approach as practiced by a financial 
institution is also effective to manage this supply -demand 
relation. In this case, it is possible to consider that 
fluctuations of the electricity price correspond to 
5 fluctuations of an interest rate in the conventional 
technique. The asset and the liability may include one 
with high uncertainty and risk or one with low uncertainty 
and risk. Therefore, it is possible to properly maintain 
the portfolio in response to the fluctuation of the 

10 electricity price (to maximize an expected return in 
response to risk tolerance) by means of appropriately 
combining short-term contracts and long-term contracts, or 
high-risk contracts and low-risk contracts. The power 
generation-procurement plan curve 103 0 in FIG. 7C 

15 conceptually illustrates this approach. In this case, the 
optimum procurement is conducted so as to correspond to 
the respective demands instead of adding the power 
generators in the order of low costs as described by the 
total power generation-procurement curve 104 0 in FIG. 7B. 

2 0 The demand and procurement are represented similarly by 
rectangular blocks in the drawings. In reality, however, 
the demand and the procurement need to be adjusted 
individually from the viewpoints of the risks and the 
profits. In other words, it is necessary to conduct the 

25 optimum procurement to meet the demand instead of simply 
operating the power generators in the order of low costs. 

Concrete methods to be used for achieving the 
optimum procurement include maturity ladder analysis, term 
gap analysis or duration gap analysis. The maturity 

30 ladder analysis is an analysis method to compile 
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maturities of the assets and liabilities which are 
currently possessed. This method is suitable for 

producing a database of the assets. The term gap analysis 
is a method of compiling the data of the maturity ladder 
5 analysis while dividing the data into each section of 
terms such as one month and thereby grasping gaps of 
matured amounts among the respective terms. This method 
is used to calculate sensitivity of profits when interest 
rates fluctuate. 

10 Although there may be a case in which an interest- 

rate risk needs to be evaluated in the power trading, 
description will be made herein on a method of an analysis 
of sensitivity to a price risk. An assumption will be 
made herein that a relation between the electricity 

15 procurement and the electricity demand are as shown in 
FIG. 8, for example. A change in profit when the 
electricity price is increased by 100 yen, for example, is 
easily calculated by use of FIG. 8. Since the demand arid 
the supply are formed of the same amounts and terms, the 

20 price volatility risk is cancelled herein. If the 

respective amounts and terms are different between the 
demand and the supply, the price volatility risk is not 
cancelled. Duration is a tool for easily expressing such 
a situation. 

25 In the financial field, price volatility of a coupon 

bond generally becomes larger as a remaining period is 
longer, assuming that the coupons, which are cash flows, 
have the same values. Therefore, the remaining period of 
the bond functions as a scale for the price volatility. 

3 0 Accordingly, a period weighted by present values of the 
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cash flows is the duration. The duration is defined by 
the following equation (11) : 

where, r is an interest rate, which is a yield to 
5 maturity of the coupon bond, and P is a price of the 
coupon bond where Ci is paid at tj.. 

Roughly speaking , in terms of the electricity price 
as well, the price volatility risk becomes larger as the 
term of the contract is longer or the price is higher. 
10 Therefore, it is possible to evaluate the price 
sensitivity by defining an amount corresponding to the 
duration. 

When the electricity generated by the power 
generator is deemed as the asset, a merit or a demerit of 

15 the asset can be evaluated by use of the rate of return. 
A rate of return of a risk asset such as a stock or a 
security is a sum of an income gain and a capital gain. 
Since there is no income gain in the case of the 
electricity, it is only necessary to consider the capital 

20 gain attributable to the fluctuations of the electricity 
price. When a future price is unknown, the rate of return 
can be deemed as a stochastic variable and a profit or a 
loss of the electricity asset can be expressed as a 
probability distribution function. Practically, the 

25 distribution of the rate of return is estimated by use of 
past data based on the assumption that observed values are 
independent and accord to the same distribution. When one 
period from a time point j-1 to a time point j is referred 
to as a period j , the rate of return in the period j is 
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expressed by the following equation (12) : 

Si -Si . 

J — (12) 

where, Sj is the electricity price at the time point 
j . Fluctuations of a fuel price are not taken into 
account in this case. A continuous compound rate of 
return in the period j is calculated by the following 
equation ( 13 ) : 



Rj = In 



(13) 



If one period is equal to one day, then the latter 
10 portion of the equation (13) is a logarithmic daily rate 
of return. 

As described above, the duration is the remaining 
period of the portfolio weighted by the present values of 
the cash flows. In the case of the electricity asset, or 
15 liability, the cash flows c(t) are deemed to be 
continuous. Accordingly, the duration d can be defined as 
the following equation (14) : 



d = 2l^r~T, (14) 



where, r is the interest rate, Ci(t) is a cash flow 
20 at time t of an i-th element in the portfolio of 
outstanding contracts, and Xi is a contracted term thereof. 
Therefore, the duration continuously changes as a time 
function in this case. Since the cash flow changes 
depending on the contents of the contract or fluctuations 
25 of a market price, the duration of the portfolio will be 
calculated by use of a Monte Carlo simulation. It is 
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necessary to find a probability distribution of the 
duration, depending on the case. 

In general, it is known to be possible to immunize a 
market risk in the portfolio by matching the duration of 
5 the portfolio with a period of investment (L. Fisher and 
R. Weil, Journal of Business, (1971) pp. 408-431). This 
is applicable to a case of a parallel shift of an 
interest -rate risk. However, it also holds true for the 
fluctuations of the electricity price approximately in 

10 terms of short-term risk evaluation. In the case of a 
long-term risk or drastic price volatility, it is 
necessary to decide an optimum portfolio combination by 
executing the Monte Carlo simulation. 

When a future is usable, the portfolio duration can 

15 be changed without modifying a spot position. The 
portfolio risk is increased when invested in an asset with 
long duration. In this case, it is possible to shorten 
the duration by selling a future. 

Moreover, in the financial field, the portfolio is 

20 regularly reorganized by use of a swap transaction between 
a short-term interest rate and a long-term interest rate 
and the like. Such techniques are also applicable to the 
risk management of the electricity portfolio. 

FIG. 9 shows a functional configuration of an 

25 electricity procurement plan producing unit 1004 as a 
second embodiment of the present invention, which can be 
employed in the power trading risk management system of 
the first embodiment. FIG. 10A and FIG. 10B is a 
flowchart of a process carried out by the electricity 

30 procurement plan producing unit 1004. Here, description 



- 37 - 



will be made on functions for producing an electricity 
procurement plan by combining electricity to be generated 
by owned power generators and electricity to be procured 
from a market, comprehensively managing the plan as a 
5 portfolio while comparing the plan with demand data, and 
thereby maximizing a profit while maintaining a risk of a 
profit raised by sale of electricity for a certain period 
within a tolerance. 

In a market database 100, various kinds of market - 

10 related data such as spot electricity price data, 
generation price data, contract term data, those are 
necessary for this system, interest rate data are stored. 

Since a power generation cost is expressed by a 
function of an amount of power generation, a power 

15 generation cost calculating unit 101 determines the power 
generation cost by use of a cost function 100A for a power 
plant, demand data 100B and fuel price data 100D (Step 
S101) . Although a conventional method can be used herein, 
it is also possible to use the method to be described 

20 later. 

A total demand curve producing unit 101A produces a 
total demand curve based on the demand data 100B (Step 
S101A) . 

A portfolio producing unit 102 firstly produces an 
25 electricity portfolio as a starting point (Step S102). 
For this portfolio, it is possible to assume that one or 
plural generators owned by a power company are employed at 
a given proportion in order to supply electricity for the 
whole demand, or it is also possible to assume that 
3 0 electricity to be generated by the owned generators and 
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electricity to be procured from the market at a proportion 
2:1. 

A power generation planning unit 103 produces the 
electricity procurement plan while comparing the obtained 
5 power generation cost with an electricity price 100C in 
the spot market (Step S103) . The procurement plan 
produced in this step is quantitatively balanced with the 
demand. However, the risk is not necessarily optimized at 
this stage. Accordingly, the portfolio risk, such as the 

10 VaR, is evaluated with a profit distribution evaluation 
unit 105 by use of volatility, which is historical 
volatility, which is calculated by using amounts demanded 
and past data for spot market prices to be provided from a 
database 104 (Steps S 105 and S106) . 

15 A risk tolerance judging unit 107 performs judgment 

of a risk tolerance by use of risk amount evaluation data 
(Step S107) . When this value is equal to or greater than 
the tolerance, a portfolio reorganizing unit 108 reduces 
the risk by reorganizing the portfolio (Step S108). When 

20 it is not possible to reduce the portfolio risk to the 
tolerance or below, reduction of the supply is also 
considered. On the contrary, when the portfolio risk is 
lower than the tolerance, a profit maximization judging 
unit 107-1 judges whether or not the profit is maximum. 

25 In a case that the profit is not maximum, an investment in 
another asset having a higher risk and a higher return is 
also considered in the portfolio reorganization unit 108 
(Step S107-1) . In this way, the profit is maximized while 
maintaining the risk within the tolerance. Moreover, a 

30 power generation-procurement plan executing unit 109 
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executes the procurement of electricity based on the 
produced power generation-procurement plan (Step S109) . 

A position review unit 111 reviews portfolio risks 
110A, positions HOB and current ALM data 100C everyday 
5 (Step Sill) and stores the result of review in a database 
112 (Step S112). An output unit outputs a risk-related 
index (Step S113) . Furthermore, when necessary, the 
position review unit 111 reorganizes the portfolio through 
transactions with a market and the like. When 

10 appropriate, the position review unit 111 performs risk 
hedging of the portfolio by use of a risk hedge executing 
unit 114 (Step S114) . The risk hedge executing unit 114 
can use instruments such as options, swaps, and futures. 

Next, as a third embodiment of the present 

15 invention, an electricity procurement plan producing unit 
1004 of the power trading risk management system will be 
described with reference to FIG. 11. Here, together with 
the input data 110A to HOC , emission right market data 
121 and weather derivative price data 122 are used in 

2 0 addition to the second embodiment shown in FIG. 9. Note 
that other functions of the third embodiment shown in FIG. 

11 are similar to those in the second embodiment shown in 
FIG. 9. 

A process carried out by the electricity procurement 
25 plan producing unit 1004 is as a flowchart shown in FIG. 

12 and FIG. 10B. A difference from the flowchart of the 
second embodiment shown in FIG. 10A and FIG. 10B is that, 
in this embodiment, a price of emission right of carbon 
dioxide 121 and a weather- related derivative 122 are to be 

30 included at the Step S102. Other steps are similar to 
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those of the second embodiment . 

By combining an emission right for carbon dioxide or 
the like with the existing portfolio, it is possible to 
maximize the profit while maintaining the risk . of 
5 fluctuations of the electricity prices within the 
tolerance. When the electricity portfolio includes power 
generation facilities such as a thermal power plant 
designed to emit greenhouse gases, and when emission 
rights for the gases are being traded in a market, the 

10 profit may be further optimized by stopping the thermal 
power plant and selling the emission right at the same 
time. The emission right in this case is not only limited 
to carbon dioxide, but is also applicable to other gases 
such as S0 2 and NO x . 

15 The weather derivative is normally used for hedging 

a weather risk. However, the demand and the price of 
electricity have a strong correlation with the weather. 
Accordingly, there is a possibility of reducing the 
portfolio risk by incorporating the weather derivative 

20 appropriately into the portfolio. 

Upon evaluation of a future risk of the electricity 
price, the relation with the demand is not ignorable. The 
future electricity price usually has daily or weekly 
periodicity, and a seasonal factor on a yearly basis is 

25 also observed. Since regular fluctuations are not risk 
factors, such regular fluctuations need to be 
appropriately processed. However, regularity of the 

electricity price is mostly unclear and it is normally 
difficult to process the regularity directly. On the 

30 contrary, regularity of the electricity demand is 
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relatively clear because the demand reflects social 
activities. Therefore, it is more convenient to evaluate 
the regularity of the price by use of the regularity of 
the demand and the relation between the demand and the 
5 price rather than directly processing the regularity of 
the electricity price. Nevertheless, the relation between 
the demand and the price is determined by a configuration 
of regional power sources, types of customers, quantity, 
composition, and characteristics of power systems. 

10 Accordingly, it is not possible to determine a universal 
relation in advance. 

FIG. 13A and FIG. 13B show a relation between the 
electricity price and the demand at the CalPX. FIG. 13A 
represents data in February, 1999 and FIG. 13B represents 

15 data in August, 1999. As apparent from these graphs, the 
relation is changed depending on the season even in the 
same region. This is attributable to the fact that types 
of operating power source facilities are different 
depending on demand levels. Moreover, in the relation 

20 between the demand and the price, a so-called quantity of 
reserved power source is important. 

In the meantime, FIG. 14A and FIG. 14B are examples 
of a relation between the demand and the price at Leipzig 
Power Exchange. The examples apparently show larger 

25 variation as compared to the case of California Power 
Exchange. Such variation needs to be taken into account 
for the risk evaluation as an uncertainty. 

As described above, the relation between the 
electricity demand and the electricity price, which varies 

3 0 depending on the country, region, and season, should be 
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appropriately adopted for the risk evaluation. 

FIG. 15 shows another functional configuration of an 
electricity procurement plan producing unit 1004 of a 
fourth embodiment of the present invention, which can be 
5 employed in the power trading system of the first 
embodiment.. The system of this embodiment is 

characterized in that the portfolio is composed by use of 
the relation between electricity demand and the 
electricity price. A demand-price correlation measuring 

10 unit 133 finds a correlation equation and a correlation 
coefficient by means of regression analysis using past 
electricity demand data 131 and past electricity price 
data 132. The demand-price correlation measuring unit 133 
outputs the correlation equation and the correlation 

15 coefficient thus obtained to a covariance matrix producing 
unit 134 and a periodicity component removing unit 135 
(Step S133) . 

The covariance matrix producing unit 134 produces a 
covariance matrix (Step S134) . Meanwhile, the periodicity 
2 0 component removing unit 13 5 removes a periodicity 
component of the electricity price, and then produces and 
evaluates statistical data 104 such as appropriate price 
volatility out of the price data after removal of the 
periodicity component therefrom (Step S135) . In the 

2 5 meantime, correlation data between the demand and the 

price can be also used for composition of the electricity 
portfolio in a portfolio producing unit 102 (Step S102) . 
Other functions of the system according to the third 
embodiment are similar to those in the second embodiment 

3 0 shown in FIG. 9. 
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The electricity price shows the daily, weekly, or 
yearly periodicity in accordance with changes in the 
demand. This periodicity is essentially different from a 
random change, but is a regular variation portion. For 
5 this reason, direct data processing may result in 
overestimation (or underestimation) of the variation 
portion, which is a risk. Accordingly, when evaluating 
the future risk by use of the past data, it is necessary 
to remove the periodicity component from the past data 
10 appropriately. There are various methods for removing the 
periodicity component. However, one typical practice is 
to decide coefficients (a-j, bj) by a least-square method on 
the assumption of a function form as shown in the 
following equation (15) : 



where, i is time or days and L is a length of the 
period, which is equal to 24 on the daily basis and or to 
seven on the weekly basis. As for the value m, a value 
around 12 is sufficient for the daily data and a value 
20 around three is sufficient for the weekly data. FIG, 14 
shows a result of removing the periodicity component out 
of the original data according to the above-described 
method . 



25 E, which is exposure, needs to be considered. The risk 
sensitivity E is defined by the following equation (16) : 



15 




(15) 



When evaluating the portfolio risk, risk sensitivity 



Ax 



(16) 



where, Ax is a change in the market data and AP is a 



- 44 - 



change in the portfolio value in that event. For example, 
if x is an interest rate, then E is the duration (of the 
conventional technique) . If x is a stock index, then E is 
a so-called beta value. Moreover, if x is an underlying 
5 asset of an option and P is an option price, then E 
corresponds to a delta factor of the option. The E in 
this case is a primary differential coefficient. However, 
it is also possible to further consider an amount 
corresponding to convexity of an interest -rate risk or a 

10 gamma factor of an option by taking consideration of a 
second differential coefficient. 

Next, as a fifth embodiment of the present 
invention, a risk evaluation method carried out in the 
power trading risk management system of the first 

15 embodiment will be described. The power trading risk 
management system according to the first embodiment 
described above is based on an assumption of using the 
normal distribution, which is generally used as the 
distribution, of the daily rates of return. However, as 

20 described above, the distribution of the daily rates of 
return is largely deviated from the normal distribution in 
the case of the electricity price. Such a fact may 
constitute a large risk. Accordingly, it is necessary to 
evaluate deviation of the distribution of the daily rates 

25 of return from the normal distribution at any time. 

FIG. 18A and FIG. 18B are graphs plotting variations 
with time of skewness a 3 and kurtosis a 4 by using the data 
of California Power Exchange (CalPX) . If the distribution 
of the daily rates of return is close to the normal 

30 . distribution, then a 3 ~0 and a 4 ~3 hold true. It is possible 
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to check the deviation from the normal distribution by- 
checking these values. Moreover, by calculating a 
stochastic amount JB ( Jarque-Bera) defined by the 
following equation (17) , for example, it is possible to 
5 evaluate a degree of deviation from the normal 
distribution : 

j M^ 2+ ^ ( * 4_3) 1 D; " 2(2) (i7) 

where, n is the number of the data. According to the 

statistic theory, it is evident that JB accords to a % 2 

10 distribution for the degree of freedom 2. Therefore, it 
is apparent that the distribution of the daily rates of 
return cannot be used as the normal distribution at 95% 
confidence level when JB is equal to or greater than 
5.991. The power trading risk management system of the 

15 fifth embodiment is characterized in that the risk 
evaluation is more strictly executed by use of the above- 
described values. Note that the stochastic data storage 
104 stores necessary data in advance. 

Next, as a sixth embodiment of the present 

20 invention, another risk evaluation method carried out in 
the power trading risk management system of the first 
embodiment will be described. In the risk evaluation to 
be carried out by the profit distribution evaluation unit 
105, the distribution obtained from the market is often 

25 used directly. However, such an approach is not accurate 
in a strict sense. The distribution of the daily rates of 
return must be obtained based on a price volatility model 
without contradictions . 

Accordingly, the system of this embodiment is 
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characterized in that a financial Boltzmann model (Y. 
Uenohara et al . , Proc . 5th JAFEE Int. Nat. Conf . , pp. 18, 
1999) is used as the price volatility models, that 
calculation is executed in accordance with the financial 
5 Boltzmann model upon evaluation of a variation component, 
and that to obtain a risk-neutral probability 
distribution, risk measurement is executed. 

The financial Boltzmann model is an expansion model 
of a diffusion model, which can evaluate a derivative 

10 security price relevant not only to the normal 
distribution but also to a price distribution of a i wider 
range . The Boltzmann model allows incorporation of a fat 
tail without damaging continuity. In this way, 

duplicatability is guaranteed and risk hedging of a 

15 derivative security is thereby facilitated. The financial 
Boltzmann equation is expressed as the following equation 
(18) : ' 

dP(Sj) +Sr dP(S,t) 



dt as 

+ ^dvd jj. 



Sv M dp(S ^° + A T (S, v)p(S, v, /i, 0 {1Q) 

Oi3 



-\dv'dn'Sp(S,v,n,t)K s (S,v,n' -* v,M)] 
= S(S-S 0 )5(t) 

where, P is a risk-neutral probability measure of an 
20 underlying asset S, t is time, v is a logarithmic rate of 
return, \i is a direction of price change, A T is a 
collision frequency, As is a scattering term representing 
a memory effect, and S 0 is S when t = 0. In the meantime, 
P(S, t) is expressed by the following equation (19): 
2 5 P(S,t)= \dvdw(S,V,/u;t) (19) 

The scattering term As is calculated on an 
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assumption of a function form as shown in the following 
formula (20) as the distribution of the logarithmic daily 
rates of return: 



v 

77 ex P 



r(v') 



T(y>) 



(20) 



5 Meanwhile, T (v) is a parameter corresponding to 

temperature, which is expressed by the following equation 
(21) : 

T(v) = T 0 (l + c 0 v + g 0 v 2 ) (21) 

Here, T 0/ c 0 , and g 0 are constants. 

10 FIG. 19A to FIG. 19D are calculation examples of the 

fluctuations of the electricity prices by use of the 
financial Boltzmann model. These graphs also include the 
corresponding distributions of the daily rates of return. 
FIG. 19A is a calculation example when the distribution of 

15 the daily rates of return is close to the normal 
distribution. FIG. 19B is a result when the distribution 
of the daily rates of return is largely different from the 
normal distribution. In the financial Boltzmann model, it 
is possible to find a risk-neutral probability density of 

2 0 a distribution in a wide range by selecting T 0 , c 0 , and g 0 
described above . 

Since the financial Boltzmann model can treat the 
distribution deviated from the normal distribution as 
described above, it is apparent that the financial 

25 Boltzmann model approximates the actual daily rates of 
return more properly than the normal distribution, and is 
therefore suitable for describing the fluctuations of the 
electricity prices. 

FIG. 20A and FIG. 20B are results of obtaining 
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accumulated distribution functions regarding the 
distribution of the daily rates of return obtained by the 
financial Boltzmann model and regarding the normal 
distribution in the case of FIG. 19A and FIG. 19B. For 
5 example, when evaluating the value at risk at 95% 
confidence level, reference should be made to the daily 
rates of return in the drawings where the accumulated 
distribution is equal to 0.05. FIG. 20B is an expanded 
view of FIG. 20A. In terms of values of the daily rates 

10 of return corresponding to the accumulated distribution at 
5%, which corresponds to the VaR, the values are not 
largely " different between the financial Boltzmann model 
and the normal distribution. However, the financial 
Boltzmann model shows a higher probability of causing a 

15 larger loss. This aspect corresponds to the actual 
situation as shown in FIG. 5. 

Although the above -described results correspond to 
the case where the VaR figures are not largely different 
by coincidence due to the difference in distribution, the 

2 0 VaR figures may be largely different when the confidence 
level of the VaR accounts for 99% or 90%. 

In the above-described results, a method of 
selecting the parameters T 0/ c 0 , and g 0 has not been 
specified in the calculation of the financial Boltzmann 

25 model. For 'example , a conceivable method is to select the 
parameters similarly to fitting to the market data. 

Alternatively, it is also possible to select T 0 , c 0 , 
and g 0 similarly to a method practiced in the financial 
and securities field. FIGs . 21A and 21B show relations 

30 between the logarithmic daily rates of return and the 
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demand, or a traded volume, in terms of the CalPX 
electricity prices and the stock prices of Company A. 
Despite ambiguity, it is still possible to observe similar 
relations from these graphs such that the demand or the 
5 traded volume is increased along with an increase in the 
rate of return. At least, clear characteristic difference 
is not recognized between the electricity price and the 
stock price. Therefore, it is deemed possible to apply a 
method based on a memory effect of the daily rates of 

10 return, which is used in the financial and securities 
fields (Y. Uenohara et al . , Proc . 5th JAFEE Int. Nat. 
Conf . , pp. 18 , 1999) . 

FIG. 22 and FIG. 23 show outlines of this method. 
FIG. 22 shows a relation between a daily rate of return on 

15 a current day and a daily rate of return on the preceding 
day. It is apparent from this graph that the shape of the 
distribution of the daily rate of return on the current 
day is more flattened as the daily rate of return on the 
preceding day is greater. The parameters T 0 , c 0/ and g 0 

20 can be determined by means of fitting with equations 18 
and 19. FIG. 23 shows a result of fitting. In the above- 
described method, the risk management with higher accuracy 
can be achieved by using the financial Boltzmann model for 
the risk evaluation of the electricity portfolio. 

2 5 Next, a power trading risk management system 

according to a seventh embodiment of the present invention 
will be described. The power trading risk management 
system of this embodiment is characterized by its profit 
optimization process carried out in a portfolio producing 

30 unit 102 on producing a portfolio as shown in FIG. 9. 
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Therefore, the functional configuration of the system of 
this embodiment is similar to that of the first 
embodiment. The portfolio producing process is also 
similar to that of the flowchart shown in FIG. 10A and 
5 FIG. 10B. 

FIG. 24 shows a functional configuration of the 
portfolio producing unit 102 of this embodiment. Here, 
the electricity generated by a plurality of power plants 
(1) to (n) with different costs are allocated to a 

10 plurality of customers (1) to (m) with different prices 
and different amounts demanded. 

When deciding an optimum power generator output of 
each power generator in response to a given electricity 
demand, power generation outputs have been heretofore 

15 allocated so as to equalize incremental fuel costs in 
accordance with the equations (9) to (13) as described 
above. However, when the power trading is liberalized, 
there is a possibility that the electricity price varies 
depending on the time or a counterpart. The allocation of 

2 0 the power generator to the demand is decided by the ALM 

approach and the methods of the portfolio optimization as 
described in the first to sixth embodiments. In this case 
as well , the amount of power generation with each 
generator needs to be determined by use of a fuel cost 
25 function. Nevertheless, whereas the generation power has 
been conventionally allocated so as to minimize the fuel 
costs, the generation power is allocated so as to maximize 
the profit in the present invention. Here, for the 
purpose of simplification, consideration will be made on 

3 0 the simplest case where there are two power generators and 
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two customers, and a power generator 1 supplies the 
electricity to a demand 1 while a power generator 2 
supplies the electricity to a demand 2. In this case, the 
equation 11 can be replaced by a question of deciding P x 
5 and P 2 for maximizing the following formula (22) under the 
condition of Pi + P 2 = Po : 

f r [/>S l (/>,0 + P 2 S 2 (P 2 j)]dt 

(22) 

where, T is a term of a contract, P x and P 2 are 
amounts of electric energy of the demand 1 and the demand 

10 2, and Si and S 2 are prices per unit electricity for the 
demand 1 and the demand 2, all of which constitute 
functions between amounts demanded and time that generally 
vary depending on customers. Meanwhile, fi and f 2 are 
power generation cost functions. 

15 Here, for the purpose of simplification, an effect 

of the present invention will be indicated based on an 
example in which Si and S 2 do not rely on demand levels and 
are constant at any time, which means at no risk. In this 
case, the formula (22) is modified to the following 

20 formula (23), because the formula (22) does not rely on 
the time and the integrals are removable: 

p x s x +p 2 s 2 -\f x m+f 2 (p 2 )] (23) 

Accordingly, assuming that 
F = P t S t +P 2 S 2 -[/i(^) + / 2 (P 2 )]+A(P 0 -/> -P 2 ) (24) 
25 the condition for maximizing the profit is expressed by 
the following equations (25) : 
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BP, dP, 

dP 2 2 dP 2 



(25) 



By resolving the equations (25) on the assumption 
that f 1 (Pi)=a 1 Pi 2 and f 2 (P2) =a 2 P 2 2 , the following equations 
(2 6 ) are obtained : 



R = 



P 2 = 



_ (S l -S 2 )/2 + a 2 P 0 
a x + a 2 
{S 2 -S t )/2 + a t P 0 



F = 



{S x -S 2 ) 2 /A + a x a 2 P 0 2 



(26) 



These results are identical to the results of the 
above-described equations (11) to (13) when Si=S 2 . The 
foregoing results reflect the case where Si and S 2 do not 
depend on the time or the demand levels. However, when S x 
10 and S 2 depend on the time or the demand levels, it is 
obvious that the results of this method are dif f erent . f rom 
the results of the conventional methods, and that the 
results of the present invention will bring a larger 
profit . 

15 Description has been made above regarding the case 

without the risk. In reality, however, S x and S 2 vary with 
time, in other words, there is a risk. In this case as 
well, optimum solutions can be found by use of the formula 
(22). When variances (a 2 ) of random fluctuations between 

20 Si and S 2 are equal to each other, the entire risk is 
decided solely by P 0 and a. Therefore, P 0 should be 
decided so as to set the risk, such as the VaR, within the 
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tolerance and then Pi and P 2 should be optimized so as to 
maximize the profit. When S x and S 2 have different risks, 
it is necessary to calculate risk values and profits in 
terms of numerous combinations of Pi and P 2 . Similarly, 
5 optimization is theoretically possible even if the number 
of demands is increased. However, computational 

complexity is drastically increased. In this case, it is 
possible to conduct calculation in reasonable time by use 
of dynamic programming method or the like. Meanwhile, fi 
10 and f 2 become time functions when fluctuations of fuel 
costs are taken into account; however, a similar method is 
applicable . 

FIG. 25 shows one example of a screen layout of a 
risk management system for an electricity portfolio 

15 according to the present invention. Reference numeral 301 
denotes targeted period, reference numeral 3 02 denotes new 
assets, reference numeral 303 denotes assets inputted in 
the past, reference numeral 3 04 denotes positions at 
present, reference numeral 305 denotes relations between 

20 prices and demands, reference numeral 306 denotes a 
distribution of a daily rates of return, reference numeral 
307 denotes an evaluation portion of VaR, reference 
numeral 308 denotes an option period, reference numeral 
309 denotes an option type, reference numeral 310 denotes 

25 fluctuations of demands in the past, and reference numeral 
311 denotes fluctuations of prices in the past, 
respectively. This system is designed to allow a user to 
select either an operation using the Black-Scholes 
equation or an operation using a Boltzmann model. 

30 The above-described power trading risk management 
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system of present invention is realized by a single 
computer, or by a network system including, a plurality of 
computers connected through a network and more computers 
dispersed in many locations which are connected through an 
5 information network. In addition, the technical scope of 
the present invention also encompasses a software program 
which is installed in a single computer or in a computer 
network system for achieving the functions of the system. 



- 55 - 



